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Abstract 


A canonical Hamiltonian formalism is derived for a class of Ermakov sys¬ 
tems specified by several different freqnency fnnctions. This class of systems 
comprises all known cases of Hamiltonian Ermakov systems and can always 
be rednced to qnadratnres. The Hamiltonian strnctnre is explored to hnd 
exact solntions for the Calogero system and for a noncentral potential with 
dynamic symmetry. Some generalizations of these systems possessing exact 
solntions are also identihed and solved. 


1 Introduction 

Ermakov systems have been intensively stndied since the late sixties both in 
view of their nice mathematical properties and of the application potential 
of their celebrated invariants. More recently, the identihcation of additional 
strnctnres in the original Ermakov-Pinney system and in its generalized ver¬ 
sion, the Ermakov-Lewis-Ray-Reid (ELRR) system (see e. g. 0] and @] and 
references therein for an npdated appraisal of the snbject and its applications) 
has called for extra attention. After the early scrntiny of the symmetry prop¬ 
erties of the ELRR system []^ and the exploration of several generalization 
schemes P, the attention has been, more recently, centered on the joint 
existence of a second independent constant of motion or a Hamiltonian for 
snb classes of Ermakov systems |^-[§. This particnlar topic is important 
per se becanse, for Ermakov systems, the existence of a second constant of 
motion nsnally implies, as shown in section 2, complete integration. Most 
importantly, however, the Hamiltonian strnctnre is fnndamental in varions 
contexts of Physics like qnantization and pertnrbation theory. 

A generic ELRR system for two independent variables is given by eqna- 
tions (B-i below. This paper considers Ermakov systems for which ut is 
generalized to be a fnnction of the variables x and y besides the time t. By 
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choosing the function / and g appropriately, it is possible to fit the equations 
into a Hamiltonian formalism and eventually derive another constant of mo¬ 
tion, distinct from the commonly known Lewis-Ray-Reid invariant (LRRI) 
given by equation (^. When two integrals exist for this Hamiltonian Er¬ 
makov system then we show that the equations are completely integrable. 

A Hamiltonian structure for the ELRR system can be enforced by spe¬ 
cializing the arbitrary functions /, g and u. For w = uj{t), that is, when uj is 
a function of t only, the Hamiltonian constraint involves only / and g. This 
specialization has already served to important applications 0 where the po¬ 
tential in the resulting Hamiltonian is essentially built with what remains of 
the functions / and g. In a more general situation where uj depends also on 
the dynamic variables, the Hamiltonian constraint results less restrictive and 
applies for a wider class of admissible systems. 

In this paper we consider the Hamiltonian property of the generalized 
ELRR system in which uj may depend not only on time but also on x and y. 
In section 2, a quadratic form in the momenta is proposed for the Hamiltonian 
and the consequences of this choice are explored analytically. In this way, 
all known cases of Hamiltonian Ermakov systems in two spatial components 
are recovered and generalized. The general class of nonlinear systems thus 
determined is, in addition, shown to be exactly integrable. This remarkable 
fact encounters application in several areas of Physics. In section 3 we apply 
the technique in two different situations which in our understanding illustrate 
its applicability both in recovering results already known in the literature and 
in identifying new exactly solvable models. Among the new results we quote a 
modified version of the Calogero potential 0 and a variation of a noncentral 
Hartmann potential, known to possess dynamic symmetry 0-0. 


The most general system in two conhguration variables that qualihes for 
an Ermakov or ELRR system is usually written as 


x + uj'^x = -^f{y/x), 
yx^ 

y + uj'^y = -^g{x/y), 
xy^ 


( 1 ) 

( 2 ) 


where the over dot stands for time derivative, / and g are arbitrary functions 
of their arguments and uj is an arbitrary function of time t and of x, y and 
their time derivatives of hrst and higher orders. For practical reasons, we 
shell consider only uj = uj{x, y, t). 
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The system of equations possesses the Lewis-Ray-Reid invariant 

0 

I =\{xy- yxf + f{s) ds + J g[s) ds . (3) 

The invariant persists for arbitrary dependence of a; on a; and y |]^, |]. 
We can therefore merge o;^ and in a single function and redefine / according 
to the following rules: 



= u}‘^{x,y,t) - ^g{x/y), 

(4) 


xy-^ 


/(s) ^ 


(5) 

9{s) ^ 

0. 

(6) 


These redefinitions simplify future considerations and cast the ELRR system 
and the LRRI in the more compact forms 

x + Vt^x = -^F{y/x), (7) 

yx^ 

y + ^^y = 0 , ( 8 ) 

I =\i.xy - yxf + F{s)ds. (9) 

The transformations (ihl) also indicate that the conventional ELRR sys¬ 
tem comprises only two arbitrary functions and not three as implied by the 
traditional notation. Of course, all previous results found in the literature, 
obtained in the standard notation, remain true. Notice, however, that we 
can always consider the Hamiltonian property of ELRR system in the form 
(0^ without any loss of generality. 

2 Hamiltonian Formalism 

As already mentioned, u in (M) can be a function of time and of any 
combinations of the dynamic variables x, y and their time derivatives of 
arbitrary order. In this work we shell consider the Hamiltonian property 
of the ELRR system for which the frequency function is allowed to depend 
not only on time (as usual) but also on x and y. In this case, the resulting 
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constraint becomes less restrictive and we can find a mnch wider class of 
ELRR systems satisfying the Hamiltonian property. 

We shall consider for the Hamiltonian of the Ermakov system (]^-^) the 
fnnction 

H = -^pI + Bp^Py + -^pI + V(x, y, t ), ( 10 ) 

where H, B and C are nnmbers snch that AC — 7 ^ 0 and V(x, y, t) is a 

potential fnnction depending on time and on the spatial variables. 

The ansatz (pRl) is jnstified in Donglas theory for two dimensional La- 
grangian systems |]^, where it is shown that the coefficients of the qnadratic 
terms in the velocities in a Lagrangian are constants of motion, at least for 
velocity-free force fields. In particnlar, these coefficients can be taken as nn- 
merical constants. As one can show, the addition of a term linear in the 
momenta does not alter the generality of the description. Finally, two cases 
of Hamiltonian Ermakov systems known in the literatnre are of this proposed 
form P, P] . 

We now impose that the canonical Hamilton eqnations generate the ELRR 
system (M- The Hamiltonian fnnction (|T^) generates the following eqna¬ 
tions for the motion of the system: 


X 

= Apx + Bpy , 

( 11 ) 

y 

= Bpx -h Cpy , 

( 12 ) 


dV 

(13) 

Px 

dx 


dV 

(14) 

Py 

dy ■ 


This first order system of eqnations in {x,y,px,Py) can be easely recast in 
the eqnivalent second order system of eqnations for (x, y), 


■■ ,dV av _ 

x + Vt^x = —A— - B ——h , 

ox ay 

OX oy 


(16) 

(16) 


where the terms proportional to were conveniently added to both sides. 
The comparison of eqnation (IT) with (H) leads to the conclnsion that the 
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admissible frequencies must satisfy 




y dx 


dy 


( 17 ) 


Also, the comparison of equation ([l5|) with equation (|^), when is given 
by (O), shows that the potential must obey the linear hrst order partial 
differential equation 


dV dV 1 

(Bx -Ay)— + (Cx -By)—= -^F(y/x). 


(18) 


The characteristic equations associated to 

dx dy 


can be written in the form 
x^dV 


Bx — Ay Cx — By F{y/x)^ 

and have - as can be easely checked - the general solution 

1 . , 1 

q 


F = + - J F(s')ds'. 


(19) 


( 20 ) 


Here A{q,t) is an arbitrary function of its arguments, and the new variables 
q and s are dehned in terms of the dynamic variables and the parameters of 
the Hamiltonian as 


q = Ay^ — 2Bxy + Cx ‘^, (21) 

s = y/x. (22) 

It is also convenient to dehne the function 

^{s) = As‘^-2Bs + C, (23) 


so that q = x‘^^{s). 

In the conventional notation, the resulting Hamiltonian ELRR system 
implied by the admissible frequencies and potentials, now reads 


dA 

X + p^x 
oq 

dA 

V + P^V 


-^fiy/x) , 

yx^ 

(24) 

-\9{.x/y), 

xy^ 

(25) 
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where 


/(«) 

^(1/s) 


2/5^ /’ - -B)F(s), 

2p 5 /’ + y(Bs - C)F(s), 


( 26 ) 

(27) 


and 


p = AC-B\ (28) 

There remains two arbitrary fnnctions in the Hamiltonian ELRR system, 
namely A(g, t) and F{y/x). In fact, / and g in are determined by the 

single homogeneons function F, and one can check directly that they satisfy 


s{Bs -C)^ + {C + 2Bs)f = {As 
as 


B)^ + (A+—)S. (29) 

as s 


It is interesting to compare this result with those in the literature. Cervero 
and Lejarreta’s Hamiltonian Ermakov systems are obtained from the 
present formalism by setting A = C = 1, B = 0, and A = oj‘^{t)q in (|^) 
and (^). The resulting functions / and g do satisfy their Hamiltonian con¬ 
straint [P, which is precisely relation (|^) specialized for the appropriated 
parameter values. This formalism has already been used to study the prop¬ 
agation of elliptic gaussian beams in nonlinear, dispersive media Q. Also 
the completely integrable class of Ermakov systems determined by Goedert 
[^] is Hamiltonian and derivable from the present formalism. In this case 
the right choice is A = G = 0, R = 1, and A = 2 w‘^{q')dq'. Needless 

to say, the functions / and g resulting from this prescription do satisfy the 
integrability condition stated in p. 

According to Liouville-Arnold theorem, 2n—dimensional Hamiltonian sys¬ 
tems possessing n independent constants of motion in involution with com¬ 
pact level surfaces are integrable by quadratures [^. For these systems, the 
motion is quasiperiodic and restricted to n—dimensional tori. The present 
class of four dimensional Hamiltonian ELRR systems will posses two inde¬ 
pendent constants of motion in involution provided that the Hamiltonian 
does not depend on time. In this case, H itself is a constant of the motion 
independent of the LRRI. So, when the function A(g, f), which is the source of 
time-dependence in the Hamiltonian, does not contain f, one can expect that 
the problem is completely integrable. In fact, the level surfaces H =const. 
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and I =const. are not compact in general, bnt the differential eqnations 
nevertheless are redncible to qnadratures. Conseqnently the Hamiltonian 
ELRR system treated here mnst be inclnded in the non trivial class of solv¬ 
able ELRR systems, to which belong, for instance, some systems analyzed 
by Govinder and Leach who considered freqnency fnnctions depending only 
on time 


R 


The Hamiltonian formalism can be nsed to solve elegantly the eqnations 
of motion. Clearly x and y are not the natnral coordinates for the problem. 
Changing from {x,y) to the new coordinates {q,s), recasts the Hamiltonian 
and the LRRl in the form 


H = 2 pqpl + ^A(g, t) + ^ (30) 

/ q 

I = yHs)pl + fF(/W, (31) 


where Pq and Ps are the momenta conjngate to q and s respectively. This 
transformation deconples the dynamics and allows to treat separately the 
snbsets {pq,q) and {ps,s) of the phase space. Moreover, when A is time- 
independent, H is also an invariant. In this sitnation, we can proceed in a 
manner similar to that used in the energy integral method of standard clas¬ 
sical mechanics and split the problem in two separable ordinary differential 
equations. 


[dq/dtY = 4:p{2qH — qA{q) — 21) , (32) 

(ds/dt)^ = 2q-\t)^\s) (^I - F{s')ds'y (33) 


Equations (|3^ - |33[) can be successively solved in terms of quadratures yielding 
a formal solution to the ELRR system. In fact, the solution of (^) requires 
the knowledge of q{t) obtained from (^). 

The structure of (p^ suggests the rescaling of the time variable according 
to 

dT{t) = dt/q{t). (34) 


Such rescaling is applicable globally only when q{t) is positive definite in time 
and r results monotonic. Under this circuntances system (^2| - ^) reads 


dq 

dr 


2 


4pq‘‘ (2qH - qA(q) - 21) , 


(35) 








(36) 



2^"(s) (/- /’-F(s')*') , 


which is a decoupled set of separable equations and, therefore, reducible to 
quadratures. This procedure specihes a general solution of the problem that 
involves four arbitrary constants, namely H, I and two constants arising 
from ( 


) and (|36D 

One interesting remark concerning the Hamiltonian Ermakov system is 
the following: changing A(g) according to A(g) i—>■ A(g) + ki or according 
to A(g) I—>• A(g) + 2 ^ 2 / 5 ', where ki and k 2 are constants, is equivalent to 
changing the values of the contants if or / in (|3^) and, therefore, will not 
change the nature of the integral to be performed. This transformation can 
be explored to identify variations of a known, exactly solvable, systems that 
are also exactly solvable. In one hand the addition of ki to A does not lead 
to a relevant variation of the problem, since it only implies a change in the 
numerical value of if. On the other hand, the addition of 2 ^ 2 / 0 ' to A does 
imply a qualitative change in the potential with no relevant change in the 
calculations. We only need to replace i in the original equations according 
to 

i ^ i + ^2 . (37) 


As a concluding remark to this section, we stress that the Hamiltonian 
ELRR has been reduced to quadratures. Whether these quadratures can 
actually be performed globally is a different question to be examined in each 
particular application. Another important follow-up remark concerns the 
nonlinear superposition law []^ associated to the ELRR systems. In general, 
this nonlinear superposition law is implicit in the sense that it cannot be 
actually applied in view of the coupling between the equations. However, 
when at least one of the equations decouples, the integration can be carried 
through and the corresponding nonlinear superposition law becomes explicit. 
For Hamiltonian ELRR systems, we arrive at equations (|3^) and (|3^) which 
constitute an explicit nonlinear superposition law. This is clear since s{t) is 
constructed using q{t) obtained from a decoupled equation, namely equation 


(0- 
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3 Sample applications 

In this section we work out some sample applications of the theory. In par¬ 
ticular we analyze the Calogero potential and a super integrable example of 
a noncentral potential. Several other potentials which represent generaliza¬ 
tions of either the harmonic or the Coulomb potentials could be treated in a 
similar way at least in what concerns the dynamics of the two-dimensional 
motion obtained by projection in an appropriate plane. The super integrable 
Hartmann potential and some generalizations or variations of it (see Big) 
certainly belong to this category. A generalized version of the coupled Fin¬ 
ney equations of interest in two-layer shallow-water wave theory |Q] can also 
be solved analitically by the formalism of this paper. 

3.1 The Calogero system as a Hamiltonian ELRR 
system 

As a first example to illustrates the analytical integration of a Hamiltonian 
ELRR system, we consider the Calogero potential and its associated system 
[^], which is a one-dimensional three body problem given by the Hamiltonian 

He = \{pI+pI + pI) + ^(( 2^1 - ^ 2 ? + {X 2 - x^f + ( 0:3 - XiY) 

I 0 

9\ 92 93 

(Xi-X2y {X2-X3y (Xs-Xi)^’ 

where a,gi,g 2 and g^ are non negative constants. A rescaling of time and 
space coordinates allows to set a = 1. Moreover, in view of the translational 
invariance of the problem, we transform to the center of mass and Jacobi 
coordinates 

R = 

X = 

y = 


-(Xi +X 2 +X 3 ) , 


V2 


(xi - X2) , 


(39) 


-^(xi -FX2 -2x3). 
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The center of mass only executes free motion, and the {x, y) dynamics is 
described by the reduced Hamiltonian 


H = ]^{pI+ pD + ^ + 


9i 


‘^92 


2a;2 


.X 




+ 


25'3 


(x + \/Zyf 


(40) 


As it can be easily checked, the reduced Hamiltonian (hOD is of Ermakov 
type. Moreover, being autonomous, it is integrable. The coefficients A, B 


and C in (p!0D become, in this case, A = C = 1 and 5 = 0. Consequently, 
g = + ?/^ = and the corresponding Ermakov potential reads 


g 1 + s^ / 4^2 4^3 \ 

2 ^ 2g (1-V3s)2 (l + v/3s)2j ■ 


(41) 


By comparing this form and equation (^) we find that for the Calogero 
system, 

A = g, (42) 


and 



1 + 5 ^ ( 4^2 4^3 \ 

2 C (l-v^s)2 (l+V3s)V' 


(43) 


Equation (|3^) can be solved analytically for A given in (^) and it can be 
verihed directly that the rescaling f i—^ r is properly dehned, that is, q{t) is 
positive dehnite in time. The resulting system of equations to be solved are 
now 



4g2(2i7g - g2 - 2/), 


(44) 


(l + ^r 2/-(1 + 5^) gr 


4^2 


+ 


4^3 


(l-y3s)2 (1 + V3s)2^ 


.(45) 


It is now easy to hnd the solution g(r'), where we make the convenient 
replacement r i—^ r' = y/^r (J is strictly positive for the Calogero system): 


?(^') 


21 


(46) 


where ci is the integration constant arising from (]4^). This integration con¬ 
stant and the integration constant C2 arising from the integration of (^) 
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below can always be expressed in terms of the initial position and of H and 

/. 

The fnnction s(r) can be evalnated in closed form for a few different valnes 
of Qi- We choose gi = g 2 = ds = g, which is perhaps the most interesting 
case yielding solntions in terms of circnlar fnnctions 

s(r') = tan l^^sin"^ sin(3(r'+ ca))^ ^ , (47) 

where Ca is the integration constant refered to before. 

The last two eqnations express the parametric orbits of the problem and 
is a general solntion involving fonr integration constants, namely H, I, ci and 
Ca- This solntion is in fnll agreement with the results of Khandekar and 
Lawande [0, but are obtained in a more systematic way. 

As mentioned before, the addition of 2 g 4 ^/q to A, where g^ is a constant, 
does not alter the nature of the analytic solutions. It only requires the change 
I —>■ I + g^ in all formulae. This proves the existence of closed form solutions 
for the modified Hamiltonian 

H = + ++ + ^ 

‘2g2 _^_25(3_^ <74 

{x — y/^nY {x + y/^uY x‘^ + y'^ 


(48) 


Therefore, in the original variables, the one-dimensional three body problem 
described by the Hamiltonian 


gi 


g2 


r2 

gs 


He = \{pI +pI+pI) + ^ ((2:1 - X2Y + {X 2 - Xa)^ + (xa - Xi)^) + 


35-4 


{X1-X2Y {x2-X‘iY (x3-Xi) 2 (Xi-Xa)^ + (x2-X3)2 - 1 - (xa-Xi) 


(49) 


is also exactly solvable. He can be viewed as an integrable modiheation of 
the Calogero system. Notice that the term in g^ is a new contribution that 
did not belong to the original system. We remark that this generalization is 
possible thanks to the fact that, in Jacobi coordinates, the Calogero system 
possesses the structure of a Hamiltonian ELRR system. 
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3.2 A Hamiltonian with dynamic symmetry 

The noncentral problem described by the Hamiltonian 


h = -{pI + pI) 


a 


+ % + 


g2x 


\/x^ + y'^ y'^ ’ 


( 50 ) 


where a, gi and g 2 are positive constants, is known to possess a dynamic 
symmetry group and is separable in parabolic and polar coordinates. 
Three dimensional extensions of H are super integrable systems |^, and 
have received attention as a noncentral force problem amenable to Feynman 
quantization 


Several other two-dimensional versions of three dimen¬ 
sional super integrable models can be reduced to Hamiltonian ELRR systems. 
We selected ( 0 ) as a good example to illustrate the technique proposed in 
section 2. 

A rescaling of time and space allows to set a = 2 without any loss of 
generality. For this problem also q = + y‘^ = r'^ and it is convenient 

to introduce polar coordinates so that s = tan^. The potential, in (r, 6*) 
coordinates, now reads 


l/ = — + 


r sin^ 6 


{91 + 92 cos 9) 


Comparison of equations (^) and (|^) shows that in this case 


A(g) = —, 

r 


(51) 


(52) 


and 

/ s=tan0 

F{s')ds' = {gi + g 2 cos 9) . (53) 

sm 9 

Let us restrict considerations to the cases where the rescaled time r is well 
dehned and monotonic. A detailed analysis shows that for positive dehnite 
LRRI invariant (sufficiently high angular momentum) the trajectories never 
cross the origin. For J > 0 the variable q = never vanishes and r is 
monotonically increasing as can be calculated directly from equation (|3^). 
We note that for I > 0, g 7^ 0 is also required for the right hand side of the 
equation (^2] ) to be positive definite, a necessary condition for the existence 
of real valued solutions. We therefore consider / > 0 and reduce the original 
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problem to the set of differential equations 




= 2r‘^{Hr^ + 2 r — /) , 


(54) 


(55) 


A direct integration of (with r i—>■ r' = V^t) yields the parametric 


orbit equation 



I 


(56) 


1 — sin(r' + Cl) 



|j (l - \/l + 4/(/-9i)/j|sm(T' + , 


(57) 


where Ci, C 2 are integration constants which can be expressed in terms of 
the initial position and of H and I. For negative values of the energy, the 
motion is bounded. When the energy is positive we hnd open trajectories that 
escape to inhnity for hnite values of t'. However in the original parameter t 
this process is regular and takes an inhnite amount of time as expected on 
physical grounds. Another interesting feature of the bounded motion is the 
fact that it does not explore the entire range 6 = 0 to 6 = 2 tt. This is evident 
from equation (|57|) , which inplies x_ < cos 9 < y;+, where 



( 68 ) 


It is ease to verihe that < 1 which corresponds to the fact that the 
trajectories never visit that sector of the plane where 9 < arccosx+. It 


is also clear that for sufficiently small values of g 2 there may exist another 


excluded sector around 9 = tt. This is the case when x_ > —1 which is 
possible if and only if 21 > g 2 and gi > g 2 - 

As in the hrst example, the addition of a term 2g^/q to A(g), where g^ 
is a new constant, does not alter the calculations. The new exactly solvable 
potential is given by 
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Here we should stress that the term proportional do gs is novel and represent 
an incorporation to the original potential that preserves the integrability 
property. Again this was possible thanks to the Hamiltonian character of 
this Ermakov system. 


4 Conclusions 

A quite general class of exactly integrable Hamiltonian ELRR system has 
been identihed and solved. The basic result comes from the fact that an 
Ermakov system being Hamiltonian (with a quadratic form in the momenta) 
is exactly solvable in terms of quadratures. Also the ultimate equations 
for those systems are decoupled, a fact that leads to practical nonlinear 
superposition laws. 

Another important feature of the formalism stems from the fact that for 
each exactly solvable model there always exists a modified version of the 
problem that is also integrable analytically. This provides a mechanism to 
spawn new integrable Ermakov systems. Two examples of physically inter¬ 
esting applications were treated in detail to illustrate the practical value of 
the method. The application scope of the method, however, is much wider 
and a detailed assessment of its reach is an open question that deserve further 
study. 

Some open questions are readily identihed. A hrst open question con¬ 
cerns the expansion of the method to higher dimensional Ermakov systems. 
From the point of view of Physics, this would be most interesting, mainly 
for three space dimensions. A second open question concerns the Hamilto¬ 
nian character of ELRR systems with velocity dependent frequency function. 
Still another interesting question concerns the perturbations theory of ELRR 
systems. This completely unexplored subject has not been touched so far, 
perhaps because of the lack of an Hamiltonian structure. The results of this 
paper opens a new prospective for such issues. 
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